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Abstract. Detailed theoretical study of the mean square ra- 
dius of extensive air shower electrons has been made in con- 
nection with further development of scaling formalism for 
electron lateral distribution function. A very simple approxi- 
mation formula, which allows joint description of all our re- 
sults obtained in wide primary energy range and for different 
observation depths is presented. The sensitivity of the mean 
square radius to variations of basic parameters of hadronic 
interaction model is discussed. 



1 Introduction 

In a series of publications [A. A. Lagutin et al. (1997, 1998, 
1999), see also R.I. Raikin et al. (this proceedings)] we de- 
scribed model-independent scaling property of lateral dis- 
tribution function (LDF) of electrons in the air showers. 
These results allow to reproduce the electron LDF up to 
r ~ (25 — 30)i? m . s ., where R m . s . is mean square radius 
of electrons. In particular, it was shown that if one use the 
mean square radius as a radial scale parameter of lateral dis- 
tribution function instead of Moliere unit, then normalized 
electron LDF becomes invariable on the primary energy and 
the age of cascade. In other words, the dependence of the 
shape of lateral distribution function on energy, observation 
depth, chemical composition of primaries and even features 
of hadronic interactions is well described by the variation of 
single scale parameter R m . s .. Unfortunately, such a sensitive 
characteristic can not be estimated well directly from exper- 
imental data so far as it demands high precision measure- 
ments of lateral distribution of electron component in wide 
radial distance range (from several meters to several thou- 
sands meters from the core). Thus the detailed theoretical 
study of i? m s ., including its sensitivity to the hadronic in- 
teraction model, is necessary for the comparison of our LDF 
with experimental data of different arrays. 
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In this paper we present our latest results on the mean 
square radius of electrons in extensive air showers (EAS). 
Detailed analysis of these results shows that there is a possi- 
bility of very simple joint description of all data obtained for 
primary energies Eq = (10 14 — 10 18 ) eV and atmospheric 
depths t — (614 — 1030) g/cm . We present such formula 
here together with initial results on the sensitivity of R m . s . to 
variations of basic parameters of hadronic interaction model. 



2 Calculation methods 

In order to simulate extensive air showers of superhigh 
energies we resort to the combination of full Monte-Carlo 
treatment of hadronic cascade with analytical expressions of 
the results of numerical calculations of electromagnetic sub- 
showers. Two different approaches were used to obtain the 
characteristics of pure electromagnetic cascade initiated by 
primary photon: the semi-analytical Monte-Carlo method 
and the method based on the numerical solution of adjoint 
cascade equations. For detailed description of these meth- 
ods see A. V. Plyasheshnikov et al. (1988); A. K. Konopelko, 
A. V. Plyasheshnikov (1997); A. A. Lagutin (1993); A. A. La- 
gutin et al. (1998). It should be noted here that the main ad- 
vantage of semi-analytical Monte-Carlo method is high ac- 
curacy in tracking very distant particles, because complete 
Monte-Carlo treatment is applied for low-energetic part of 
cascade with threshold energy E t h = 0.1 MeV. On the other 
hand, the numerical solution of adjoint cascade equations 
taking into account the deflection of photons in the multi- 
ple Compton scattering is preferable for independent calcu- 
lations of moments of lateral distribution function (such as 
mean square radius) and gives very high performance in fluc- 
tuation problem and sensitivity problem. 

For hadronic part of the cascade we used simplified al- 
gorithmic generator (A. A. Lagutin et al., 1997, 1998) based 
on the quark-gluon string model (A. B. Kaidalov et al., 
1986; Yu.M. Shabelski, 1986) with extrapolation to 10 18 eV. 
Though our generator is very simple, especially in compari- 



2 



R. I. Raikin et al.: Mean Square Radius of EAS Electrons 



son with modern high-developed codes, the main features of 
interactions, that make a major impact on shower develop- 
ment (including fluctuations of full and partial inelasticity co- 
efficients) are well reproduced by our program [see A. A. La- 
gutin et al. (1998); R. I. Raikin (2000)]. Besides this code is 
well adapted for the investigation of sensitivity of results to 
variations of basic parameters of hadron-nucleus interaction 
model. 

Due to complete Monte-Carlo treatment of hadronic cas- 
cade, the computation time depends strongly on primary en- 
ergy. Nevertheless, the simplicity of hadronic generator to- 
gether with rejection of tracking muons in combination with 
analytical expressions used for electromagnetic component 
allowed us to perform very speeding calculations, that ensure 
maximum 1.5% relative statistical error in i? m . s . estimation. 

Thus we calculated cascade curves and lateral distributions 
of electrons in (10 14 — 10 18 ) eV proton-induced EAS at six 
atmospheric depths from 614 to 1030 g/cm 2 in standard at- 
mosphere. Since we implemented analytical expressions of 
normalized average LDF and mean square radius for each 
partial electromagnetic subshower taking into account fluc- 
tuations of shower size caused by both hadronic and electro- 
magnetic parts of a cascade, the mean square radius of EAS 
was calculated in the following way: 



Rm.s. — 



E (NeR 



2 

m.s. 



1/2 



(1) 



were index i denotes characteristics obtained for r h indi- 
vidual extensive air shower, the sum goes over all simulated 
EAS of fixed primary energy. 



3 Results 

The results of our calculations of the mean square radius 
of electrons in vertical extensive air showers generated by 
primary protons of E = (10 14 — 10 18 ) eV for atmospheric 
depths t — (614 - 1030) g/cm 2 are presented in Fig. 1. One 
can see from this figure, that at fixed atmospheric depth i? m . s . 
decreases with primary energy. That is, our results indicate 
narrowing lateral distribution of electrons with primary en- 
ergy in whole considered energy range. Though narrowing 
rate NR = 9R m . s ./91ogE becomes smaller when energy 
increases, it is still essential at 10 18 eV. It can be understood 
as a reaction to the decrease of ages of powerful partial elec- 
tromagnetic subshowers, which give main contribution to the 
EAS electron yield. 

The behaviour of mean square radius with atmospheric 
depth is, however, more complicated. The basic tendency is 
obvious expansion of shower in process of development re- 
alized in increasing mean square radius. Nevertheless, there 
is distinct maximum in R m . s . (t) obtained for relatively low 
primary energies, that shifts deeper into the atmosphere when 
energy increases. The analysis of our results shows that such 
a behaviour can not be explained by only the influence of 
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Fig. 1. Mean square radius of electrons in vertical proton-initiated 
EAS. Symbols - our calculational data for E = (10 14 - 10 18 ) eV 
and t — (614— 1030) g/cm 2 . Curves - polynomial approximation 



the increase of air density. There is another one mechanism 
of the narrowing of average lateral distribution related to the 
correlation between shower size [(N e )i] and width [(R m .s.)i] 
in individual showers at fixed primary energy. 

In Fig. 2 we show this correlation by plotting (R m . s .)i 
vs. (N e )i at sea level (a) and at 614 g/cm 2 (b) for indi- 
vidual showers of 10 14 , 10 16 and 10 18 eV (500 showers of 
each energy were included in data set). It is seen from part 
a) of a figure that at sea level (Rm. s .)i and (iV e )j demon- 
strate strong anticorrelation. Since relatively large and si- 
multaneously narrow individual showers of considered en- 
ergy give more essential contribution to average LDF [see 
also eq. (1)], the expansion of the average shower is slowing 
down with increasing fluctuations when atmospheric depth 
increases. As analysis shows, such anticorrelation exists 
at observation depths located deeper than level t opt , where 
fluctuations of (N e )i are minimal (according to our results, 
iopt = imax + 55g/cm 2 , where t max is depth of maximum of 
average cascade curve). Thus we should expect existence of 
such a maximum at higher energies in case of artificial show- 
ers, which continue their development below the sea level as- 
suming the further exponential growth of air density together 
with increasing fluctuations. On the other hand, for atmo- 
spheric depths above t opt such as, for example, 614 g/cm 2 
for Eq = 10 18 eV (see Fig. 2b) mean square radius correlate 
with shower size and this mechanism works conversely. 

In this connection it is necessary to note, that R m . s . esti- 
mated as (1) describes the shape of average lateral distribu- 
tion of electrons and deviates distinctly from average mean 
square radius (i? m . s .), which one can calculate using data 
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Fig. 2. Correlation between (N e )i and (i? m .s.)i in individual ex- 
tensive air showers of E Q = 10 14 , 10 16 and 10 18 eV at t = 
1030 g/cm 2 (a) and t = 614 g/cm 2 (b) 



about individual showers as 
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The difference becomes essential when observation level is 
located far from t opt . For example, at sea level 10 14 eV in- 
dividual showers are (on average) significantly wider than 
average shower ((i2 m . s .) is 15% larger than Rm.s.)- 

It is important that, though mean square radius is affected 
by number of different factors, a very simple approxima- 
tion of all our results based only on the information about 
longitudinal shower development is possible. On Fig. 3 the 
whole set of our calculational data is presented in the form 
Fr{ 3 ') = Rm.s.{s') x p /p(t), where p(t) is air density at 
depth t, po = 1.225 g/cm 3 , s' = t/(*max + 100 g/cm 2 ). 
It is seen that functions F R (s') for different primary ener- 
gies are precisely superposed with each other in overlapping 
s' intervals. The fit of Fr(s') with maximum relative error 
of 1% was obtained using method of maximum likelihood as 
Fr{s') = a + b ■ arctg(s' — 1), where a and b - free pa- 
rameters. Thus we can propose an analytical expression for 



^m.s.(-^O) t) as follows: 
Rm. s .{Eo,t) = -^7 173.0 x 



2 / t 

0.546 + - arctg — — ^ 
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m. (2) 



Our results on t meLX for proton-induced showers were approx- 
imated as 



W^o) = 740 + 65 lg(£ /10 ls eV), g/cm 2 . 



(3) 



Assuming the validity of the superposition model for 
showers initiated by a nucleus of energy Eq and mass A, we 
obtain the relation 

R m.sX E 0,t) = R ms (E /A,t), 

where upper indexes denote the type of primary particle. 



4 Discussion 

Expression (2) together with scaling parameterization of 
LDF (A. A. Lagutin, R. I. Raikin, 2001 ; R. I. Raikin et al., this 
proceedings) allows one to obtain the reliable data on nor- 
malized average lateral distribution function of EAS elec- 
trons up to r ~ (25 - 30)i? m . s . (~ 2000 - 3000 m) in wide 
primary energy interval and for any practically important at- 
mospheric depth from mountain level to sea level. Neverthe- 
less, there is always vital question related with applicability 
of any theoretical result in the field of EAS research: how 
strong it depends on hadronic interaction model used in cal- 
culations? 

In our papers [see, for example, A. A. Lagutin et al. (1997, 
1998); R. I. Raikin (2000)] it was shown that scaling property 
of electron LDF is very poorly sensitive to variation of basic 
parameters of model. Thus mean square radius carries the 
whole information about the influence of features of hadronic 
interactions on the shape of LDF as well as about mass com- 
position and shower development. Since the key parameters 
of hadronic interactions most influential for the longitudinal 
development of EAS are inelastic cross section of p-air col- 
lisions and inelasticity (J. Knapp, 1999), we have performed 
the basic test of sensitivity of expression (2) by the change 
of these parameters in our code. We used p-air cross sections 
from different well-known hadronic models (J. Knapp, 1999) 
(including MOCCA'92 cross section, which demonstrates 
highest increasing rate with primary energy), and indepen- 
dently varied inclusive energy spectra of secondaries in such 
a way that spectrum-weighted moments of inclusive spec- 
tra change by ±20% [see G.Battistoni (1999)]. Our results 
show, that within the limits defined by widespread hadronic 
models formula (2) remains valid. 

It is also important, that if one extrapolate expression (2) 
to ultrahigh primary energies (E > 10 18 eV), then values of 
variable s' corresponding to t « (920 — 1030) g/cm 2 will 
be kept within interval well presented in our data set by the 
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Fig. 3. Function Fr(s') = Rm. B .( s ') x Po/p(t), where s' = 
i/(imax + 100 g/cm 2 ). Symbols are same as in Fig. 1. Curve - 
fitting function (2) 

values calculated for moderate energies and higher observa- 
tion levels [for example, s' — 0.95 when E = 10 20 eV 
and t = 920 g/cm 2 , assuming the validity of formula (3)]. 
Though such extrapolation is rather formal, we expect that 
it can give reliable predictions at least in the energy range, 
where Landau-Pomeranchuk-Migdal and magnetic field ef- 
fects of hadron-initiated showers and also 7r°-air inelastic 
interactions are negligible. Of course, special calculations 
are needed to confirm the reliability of extrapolation of for- 
mula (2). 



5 Conclusion 

This paper represents the further development of scal- 
ing formalism (A.A.Lagutin et al., 1997, 1998, 1999) for 



the description of lateral distribution function of EAS elec- 
trons. The calculations of the mean square radius of elec- 
trons in extensive air showers have been made for primary 
energies E = (10 14 - 10 18 ) eV and atmospheric depths 
t = (614 — 1030) g/cm 2 . Very simple approximation for- 
mula for R m , s ,(E ,t), that gives the one-valued relation be- 
tween depth of maximum of average cascade curve and the 
shape of average lateral distribution of electrons is presented. 
It can be practically implemented during comparisons of ex- 
perimental data of large ground-based air shower arrays with 
results obtained by fluorescent detectors. 
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